Exact-Exchange Spin-Current Density-Functional Theory 
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A spin-current density-functional theory (SCDFT) is introduced, which takes into account cur- 
rents of the spin-density and thus currents of the magnetization in addition to the electron density, 
the non-coUinear spin-density, and the density current, which are considered in standard current 
spin-density-functional theory. An exact exchange Kohn-Sham formalism based on SCDFT is pre- 
sented, which represents a general framework for the treatment of magnetic and spin properties. As 
illustration an oxygen atom in a magnetic field is treated with the new approach. 



Density-functional method^, in particular Kohn- 
Sham (KS) methods, have developed into the most 
widely used approach to investigate electronic structures 
of molecules, clusters, surfaces, or solids. However, mag- 
netic properties and generally properties and quantities 
related to currents cannot be described with the available 
KS methods relying on standard spin-density-functional 
theory, which is based on the spin-density alone. In or- 
der to treat magnetic effects due to both spin and cur- 
rents within a density-functional theory (DFT) frame- 
work, Vignale and Rasolt therefore introduced a current 
spin-density-functional formalism almost 20 years agoSi^. 
Despite the fact that magnetic effects are ubiquitous 
and of paramount importance, see, e.g., magnetooptics, 
spintronics, or nuclear magnetic resonance spectroscopy, 
current spin-density-functional theory (CSDFT) did not 
play a significant role in practice so far. The reason is 
that sufficiently accurate approximations for the required 
but not exactly known current spin-density functionals 
for exchange and correlation are not available. 

In recent years exact-exchange (EXX) KS methods 
have been developed, which treat exchange energy and 
potential exactly and therefore do not require any ap- 
proximate exchange functionals^. These EXX KS meth- 
ods yield band structures for semiconductors and gen- 
erally one-particle spectra of electronic systems, which 
are strongly improved compared to methods based on 
the LDA or generalized gradient approximations^. A 
development of EXX current density-functional methods 
therefore seems highly promising^. 

CSDFT depends on seven variables: the electron den- 
sity, the X-, y-, z-components of the spin-density, and 
the X-, y-, z-components of the paramagnetic density 
currenliSi^. However, as demonstrated below, in general, 
also currents of the spin-density, i.e., currents of the mag- 
netization, may occur in the presence of magnetic fields. 
For a complete description of magnetic effects it there- 
fore seems desirable to go beyond CSDFT and to include 
also the x-, y-, z-components of the currents of the x-, 
y-, z-components of the spin-density, i.e., to include 9 
more variables leading to a total of 16 variables. The 
resulting formalism shall be called spin-current density- 
functional theory (SCDFT) as opposed to current spin- 
density-functional theory in order to emphasize that it 
depends on spin-current densities (The term spin-current 
density from now on shall include besides actual spin- 



current densities also standard current densities as well as 
non-coUinear spin-densities including the standard elec- 
tron density) . In this Letter we firstly introduce the basic 
SCDFT, then present an exact exchange KS formalism 
based on SCDFT, and finally consider as an illustrative 
example an oxygen atom in a magnetic field. 

The Hamiltonian operator of an electronic system in a 
magnetic field is given by 



= T + V,, 
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In Eq. Q], T and Vee are the operators of the kinetic 
energy and the electron-electron repulsion, Vext is the ex- 
ternal electrostatic potential, usually the potential of the 
nuclei, B designates a magnetic field with accompanying 
vector potential A. By the position of the i-ih electron 
is denoted, by Pi the corresponding canonical momentum 
operator. The vector cr contains the Pauli spin matrices. 
The sum in the first line of Eq. ^ runs over all N elec- 
trons. The vector S has four components, Eq being a 2x2 
unit matrix and Si, E2, and S3 being the Pauli spin ma- 
trices (Txi f^y, and (Tz- The four components of the vector 

J(r) are the density operator J2iLi ^('"^ ^ ^i) — ^o(r) 
and the x-, y, z- components of the current operator 

(5) J2zLi Px.t <5(r - Tj) + '5(r - Ti) px^, = Ji(r) etc. 
{Px,i = —id/dxi denotes the momentum operator of the 
«-th electron in x-direction.) The 4x4 matrix V(r) is 
composed of matrix elements V^i,(r) with /i, ^' = 0, 1, 2, 3 
and is given by 



fvextir) + ^ Ax{r) Ayir) ^,(r)\ 



V(r) 



V 



Bx(r) 
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(2) 



The total energy E of the system is given by 



E = T + Vee + J2 ^^''W Po-M-^) 

= T + Vee + Jdr V^(r)po(r). (3) 
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In Eq. P] , T and Vee denote the kinetic and the electron- 
electron repulsion energy, respectively, and po is the 
ground state spin-current density consisting of 16 com- 
ponents po.A"/ with — 0, 1, 2, 3. The component po.oo 
is the regular ground state electron density, po./iO with 
fi = 1,2,3 represent the x-, y-, z-components of the 
spin-density, pofiu with i' = 1,2,3 represent the x-, y-, 
z-components of the paramagnetic current of the elec- 
tron density, while po.p.v with fXjV = 1,2,3 represent the 
components of the paramagnetic currents of the spin- 
density, i.e., currents of the magnetization. The corre- 
sponding physical currents consisting of the contributions 
Po,/xi'{r)+po,Mo(r)^i/(r) describe, via the continuity equa- 
tion, changes of the magnetization in time and thus are 
meaningful physical quantities. The spin-current density 
Pq can be considered as a 4x4 matrix or, like in the sec- 
ond line of Eq. as a sixteen-component vector with 
superindex jiv. Similarly, V(r) can be treated as 4x4 
matrix, like in the second line of Eq. PP, or as sixteen- 
component vector, like in the second line of Eq. [3]. 

The reason to introduce the matrix V(r) and to rewrite 
the Hamiltonian operator according to the second line 
of Eq. pP is that this leads to the expression [5] for 
the total energy E. Starting from this expression for 
E it is straightforward to generalize the definitions and 
derivations of standard DFT. A Hohenberg-Kohn func- 
tional F[p\ can be defined via a generalized constrained- 
searcbSi^ by 

F[p] = Min (^|f + Kel*). (4) 

The minimization 0] runs over all wave functions yield- 
ing the spin-current density p, i.e., the wave functions 
are constrained to not only yield a given reference elec- 
tron density, like in standard DFT, but to yield the 16 
components p^i, of the spin-current density. The mini- 
mizing wave function 'I'[p] of is a functional of the 
spin-current density p. The minimizing wave function 
^'[po] for the ground state spin-current density po is the 
ground state wave fimction 'J'q of the system, ^I^o = '5 [po] , 
i.e., in case of the presence of magnetic fields, the ground 
state current-spin density determines via the constrained- 
search (0} the ground state wave function and thus all 
ground state properties of the system. The latter state- 
ment represents a Hohenberg-Kohn theorem for SCDFT. 

A Kohn-Sham wave function $[p] corresponding to a 
wave function ^'[p] is defined as the minimizing wave 
function in the constrained search minimization 

r,[p] - Min ($|f 1$) (5) 

yielding the functional Ts [p] of the noninteracting kinetic 
energy. As usually, exchange and correlation energy are 
defined as E^[p] = mp]\Vee\^[p]) - U[poo] and Ec[p] = 

(*[p]|r -t- VeeMp]) - ($[p]|T + VeeMp]) with the 

Coulomb energy C/[poo] given by the standard expression 
i Jdrdr' depending only on poo ■ 



The many-electron KS equation corresponding to a 
Schrodinger equation for real electrons with Hamiltonian 
operator ^ is determined by the KS Hamiltonian oper- 
ator 

Hs = f + Jdr V,(r) J(r) (6) 

with the 4x4 matrix representing the KS potential, 
given by V^(r) = V(r) -f U(r) + V^c(r). The Coulomb 
potential U and the exchange-correlation potential 
can be represented by 4x4 matrices; within the matrix 
U only the component Uoo(r) = Jdr' ^^^^ is different 
from zero and is the standard Coulomb potential. By 
straightforward generalization of the arguments of the 
standard KS formalism^, i.e., by invoking the variational 
principle for the KS system and the real electron system 
and by comparing the resulting Euler equationai^, it can 
be deduced that the components 'Vxc^i^ of the exchange- 
correlation potential "Vxc are given by 

Vx..Mr) = • (7) 

In the functional derivatives Q all 16 components of 
are treated as independent. As a result all 16 components 
Vxcfiv and subsequently Vs^fj,,y, including components (or 
combinations thereof) corresponding to gauge transfor- 
mations, are determined within the presented formalism. 
The gauge of all quantities is determined by the gauge 
of the external vector potential A. (A modified version 
of the formalism chooses the Coulomb gauge throughout 
and treats only the transversal components of the spin- 
currentsS.) While the KS wave function <&[po], corre- 
sponding to the ground state V&o = ^' [po] of a real system, 
is uniquely defined by the constrained-search |SJ), the 
functional derivatives Q uniquely define the exchange- 
correlation potential Yxc and subsequently yield the ef- 
fective KS potential Yg if they are evaluated at po. Like 
in the standard KS formalism, it is assumed that the 
functional derivatives Q exist and that the spin-current 
density po is noninteracting w-representable. 

The components Vxc,^i/ and subsequently Vs^^i, with 
fijV — 1,2,3, i.e., the components coupling to actual 
spin-currents, in general, differ from zero despite the 
fact that the corresponding components in the exter- 
nal potential V of Eq. Q are zero. A similar situa- 
tion occurs in the standard spin-density formalism^. If 
the spin-density of the system is spin-polarized then the 
KS exchange-correlation potential and the effective KS 
potential are spin-polarized although the external po- 
tential is not spin-polarized. Because the external po- 
tential does not contain a contribution coupling to the 
spin-polarization it is, in principle, possible to treat spin- 
polarized systems within the basic (non-spin-polarized) 
KS formalism. However, for approximate exchange- 
correlation functionals and thus in all practical applica- 
tions, a spin-polarized treatment is more preferable for 
spin-polarized systems. Similarly it is possible in prin- 
ciple, to treat within the original current spin-density 



3 



formalism systems, which are subject to magnetic fields 
because the external potential does not couple to spin- 
currents, i.e., because V^i, = for ^, = 1,2,3. How- 
ever, in the presence of magnetic fields, in general, spin- 
currents occur (see below) and therefore a treatment 
within the spin-current density formalism presented here 
seems to be preferable. 

If spin-orbit interactions are considered then the com- 
ponents V^i, for fi,!/ = 1,2,3 of the external potential 
no longer are zero and the original current spin-density 
formalism is no longer applicable and it is required to 
switch to the spin-current density formalism presented 
here. Indeed a generalization of the presented formalism 
that enables the treatment of spin-orbit interactions is 
possible and will be published elsewhereifi. 

In order to apply the presented SCDFT, suitable 
exchange-correlation functionals are required. To that 
end, we concentrate on the exchange functionals and 
treat those exactly. Like in the standard KS formalism, 
the evaluation of the exchange energy is well-known from 
Hartree-Fock methods and straightforward. By gener- 
alizing the treatment of the standard exact-exchange 
KS approachii also the exact exchange potential Va; of 
SCDFT becomes accessible. The chain rule of functional 
differentiation yields 



dr' 



__ SMr") 
ba{r") SVsMr') 



(8) 

In Eq. (jHJ, 4>a stands for the KS orbitals, complex- valued 
two-dimensional spinors, building the KS determinant. 
The right hand side (r.h.s) of Eq. © due to the sum 
over kX contains 16 terms. Furthermore, for each of the 
16 components of V^: an equation ^ arises. These can 
be combined in the 16x16 matrix equation Va;(r) = 
/ dr' Xj^(r,r') t(r'), which is turned into 



dr' X,(r,r') V,(r') = t(r) 



(9) 



Here, Xs is the spin-current response function of the KS 
system, a 16x16 matrix with matrix elements 



occ.unocc 



{(l)a\<y,,X{r)\(t>s){(t>s\(^KJ\{r')\(f)a) 



(10) 



The components t^i^ of the 16-dimensional vector t are 
given by the terms in the square brackets on the right 
hand side of Eq. ©. The functional derivative ^^^"^g,,^ 
contained in the tfj,i, can be evaluated straightforwardly 
because the exchange energy is known in terms of the KS 
orbitals, see, e.g., Ref.Q. Eq. @ represents the EXX 
equation of SCDFT; it is an integral equation, which 
needs to be solved in each iteration of the SCDFT KS 
self-consistency process in order to obtain the exchange 
potential V^,. 



We implemented the presented EXX-SCDFT approach 
within a plane-wave solid state code. The quantities 
fj.u,K\j Vx^fivi and tfj^u are expanded into an auxiliary 
plane wave basis set. The ^s,pi/,k;A then turn into ma- 
trices of the dimension of the auxiliary basis set, while 
the Vx.p.f and turn into vectors of the same dimen- 
sion. The complete response function then turns into a 
matrix consisting of 16x16 submatrices ^s./xi/.kA and the 
complete exchange potential and r.h.s of Eq. ^ become 
vectors consisting of 16 subvectors Vx,fi,u and tfj_^. The 
integral equation (PJ turns into a matrix equation and 
can easily be solved. 

As an illustration we consider an oxygen atom subject 
to a constant magnetic field directed along the z-axis. To 
that end, we put the oxygen atom in the centers of simple 
cubic supercells of 9, 10, and 11 a.u. length. Plane wave 
cutoffs of 14, 20, and 24 Rydberg for the orbitals and of 
4, 6, and 7.2 Rydberg for the response function, the ex- 
change potential, and the r.h.s of Eq. lO are chosen. All 
results shown below are converged with respect to super- 
cell size and plane wave cutoffs and are invariant with re- 
spect to gauge transformations of the external vector po- 
tential A, within the numerical accuracy determined by 
the plane wave cutoffs. The oxygen Is electrons are taken 
into account via separable EXX pseudopotentialsii. 

In the absence of magnetic fields and spin-orbit interac- 
tions the electronic ground state of oxygen is a nine-fold 
degenerate '^P-state of valence configuration (2s^,2p''). 
The spin and spatial angular momenta can be coupled to 
give ^Po, '^Pi, and ^P2 LS-coupled states. In the pres- 
ence of a constant magnetic field in z-direction, the ^P2- 
state with total magnetic quantum number Mj=-2 is 
the ground state. This state is a one-determinantal state 
with occupied p-orbitals p_i _i/2, Po,-i/2i P+1.-1/21 and 
p_i _i_i/2- (First and second subscripts, mi and mj, de- 
note magnetic quantum numbers for spatial and spin an- 
gular momentum, respectively.) A KS treatment of this 
state even in the limit of zero magnetic field leads to 
a symmetry breaking KS wave function and a lifting of 
the degeneracy of the p-orbitals because of the partially 
filled p-shell accompanied by a nonspherical but cylindri- 
cal electron density^ii^. Correspondingly, the symmetry 
of the KS Hamiltonian operator is reduced from spher- 
ical to cylindrical symmetry around the z-axis. In the 
presence of the magnetic field the latter is the actual 
symmetry of the system. 

In the absence of spin-orbit interactions the change AE 
of orbital energies (and, below, of the total energy) due 
to the magnetic field of strength B is described according 
to the Paschen-Back formula AE — (3 B [mi 2ms) 
not according to the Zeeman formula. Here /3 stands 
for the Bohr magneton. The values of AE/{f3 B) for 
the curves in Fig. ^ have to equal and, indeed, equal 
[mi + 2ms)- For the special case of the '^P2-state with 
Mj — —2 the Paschen-Back and the Zeeman formula 
AE = (3 gj B Mj yield the same energy splitting with 
respect to the magnetic field strength. It is thus possible 
to obtain from our calculation the Lande factor gj of 
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Magnetic Field [T] 



FIG. 1: Eigenvalues of oxygen KS orbitals Sme.ms andp, 
as functions of magnetic field strength 



In Fig. 121 the electron density poo and its current poi 
in x-direction are displayed along the y-axis. As usual 
in pseudopotential calculations poo exhibits a minimum 
at the nucleus. The displayed poi corresponds to a cur- 
rent around the z-axis. The currents po2 and po3 are 
zero. Thus the magnetic field correctly aligns the spatial 
angular momentum along the z-axis. In Fig. 13 the spin- 
density P30, i.e., the magnetization along the z-direction, 
and its current p3i in a::-direction are displayed along the 
y-axis. The spin densities pio and p2o are zero, thus also 
the magnetization is correctly aligned along the z-axis 
by the magnetic field. Most important, Fig. O clearly 
shows that the spin-current pai differs from zero. This 
demonstrates that spin-currents, which can be correctly 
described by the presented SCDFT, are far from negligi- 
ble in real systems. 



oxygen despite the neglect of spin-orbit interactions. The 
resulting value of 1.5 agrees with the experimental one. 
This demonstrates that the presented SCDFT approach 
works correctly in practice in the sense that it correctly 
generates spin and angular momenta and aligns them 
correctly with the magnetic field. (Note that in our plane 
wave supercell method spin and angular momenta as well 
as their coupling and alignment are not enforced by the 
basis set and by occupation numbers but are an outcome 
of the SCDFT treatment.) 
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FIG. 2: Electron density poo of oxygen and its current poi 
in a;-direction displayed along the y-axis for a magnetic field 
strength of 10000 T 
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FIG. 3: Spin density pao of oxygen and its current pai in 
x-direction displayed along the y-axis for a magnetic field 
strength of 10000 T 



In summary, we believe that the presented formalism 
and method provide a general and sound basis for the 
treatment of all sorts of magnetic and (non-colinear) spin 
effects in a density-functional framework. Furthermore, 
as mentioned above, also spin-orbit interactions can be 
included in a natural way in the formalism.^". 
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